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ABSTRACT 



, Received ?; accepted ' 

Q 

Context. Radial velocity (RV) measurements from near-infrared spectra have become a potentially powerful tool to search for planets 
around cool stars and sub-stellar objects. As part of a large survey to characterize M-dwarfs using NIRSPEC at Keck II, we obtained 
' spectra of eight late M-dwarfs (spectral types M5.0-M8.0) during two or more observing epochs per target. These spectra were taken 

I with intermediate spectral resolving powers (R ~ 20, 000) in the y-band. 

• , Aims. We search for relative radial velocity variability in these late M-dwarfs and test the NIRSPEC capability of detecting short period 

i-C ■ brown dwarf and massive planetary companions around low-mass stars in the y-band (» 1.25 /jm). Additionally, we reanalyzed the 

data of the M8-type star vBlO (one of our targets) presented in Zapatero Osorio et al. (2009), which were obtained with the same 
JL \ instrumentation as our data. 

t*^ , Methods. To achieve a precise RV measurement stability, the NIRSPEC spectra are self-calibrated by making use of the telluric 

■ absorption lines, which are present in the observed spectra and used as a long-term stable reference. In the modeling process a multi- 

^ ' parameter -optimization is employed to generate an accurate description of the observation. The telluric lines allow us to model the 

instrumental profile of the spectrograph and the determination of the Doppler shift of the stellar absorption lines. 
Results. For the entire M-dwarf sample, we do not find any evidence of relative RV variations induced by a short period brown dwarf 
or massive planetary companion. The typical RV precision of the measurements is between 180 and 300 m s"', which is sufficient to 
^ ' detect hot Neptunes around M-dwarfs. Also, we find that the spurious RV shift in Zapatero et al. (2009) of the star VB 10 was caused 

^Nj ' by asymmetries in the instrumental profile between different observing epochs, which were not taken into account in their analysis. 
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1. Introduction the MO to M9.5 spectral types. They exhibit prominent absorp- 

__l tion features corresponding to strong neutral atoms, H2O, FeH, 

^ The search for extrasolar planets has led to more than 700 con- ^q, and TiO. Owing to the low masses of these objects, 

. . , firmed discoverieil by using all detection techniques. Up to j^e reflex motion of the host star due to the gravitational pull of 

> . now, most of them have been detected by means of the ra- ^^e extrasolar planet is higher and more easily detectable than 

' dial velocity (RV) technique using high-resolution spectrographs ^^^^ ^^^^^^^ host stars. Since M-dwarfs are very cool stars 

^ (7? = A/AA > 40, 000) at optical wavelengths. Most discoveries j„ comparison with solar-type stars, short period planets would 

^ are giant gaseous planets (typically hot Neptunes and Jupiters) ^^^.^ jjj^gjy situated in the habitable zone. 
■ " " ' of short periods (of a few days) around stars of spectral types F, 

G and K. M-dwarfs emit most of their energy around 1.1 - 1.3 //m. 
As potential hosts to rocky planetary companions, M-dwarfs in the near-infrared (NIR), while they appear very faint at opti- 
have become increasingly popular as targets for RV searches cal wavelengths. First attempts to measure RV variations among 
(e.g. Endl et al. 2006, Charbonneau et al. 2009, Mayor et al. very cool M-dwarfs at NIR wavelengths were done by Martin 
2009, Zechmeisteret al. 2009). Very cool stars such as M-dwarfs et al. (2006). They achieved a RV precision of around 300 
are the most abundant type (~ 70%) of stars in the solar neigh- m s"' for the M9.5-dwarf LP944-20 by using the spectrograph 
borhood and the Milky Way in general (Henry at al. 1997). The NIRSPEC, which is mounted on the Keck II telescope in Hawaii 
effective temperatures and masses of M-dwarfs, respectively, are (McLean et al. 1998). Recently, several research groups have re- 
in the range 3700 to 2200 K and 0.5 to 0.07 solar masses for ported high-precision RV measurements taken in the NIR with 

CRIRES (Kaufl et al. 2004), mounted at the UTl/VLT in the 

Send offprint requests to: frodler@iac.es Paranal Observatory of ESO in Chile. Bean et al. (2010a) con- 

' The Extrasolar Planets Encyclopedia; |http://www.exoplanet.eu] ducted high-resolution spectroscopic data of over 60 M-dwarfs 

2OII-N0V-I5 (spectral types M4-M9) and used a NH3 gas cell spectrum as 
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a stable reference, and report an RV precision of better than 
5 m s"^ Figueira et al. (2010) took observations of the plane- 
tary candidate TW Hya and achieved a RV precision better than 

10 m s ' by adopting telluric lines as a stable reference. Blake 
et al. (2010) report RV measurements of 59 M- and L-dwarfs 
using the Keck/NIRSPEC spectrograph, with the aim to detect 
low-mass companions. They made use of strong CO absorp- 
tion features around 2.3 fim in M- and L-dwarfs and achieved 
RV precisions between 50 and 200 m s Tanner et al. (2010) 
report preliminary results of a late M-dwarf survey by using 
Keck/NIRSPEC with RV precisions between 150-300 m s 

In 2009, Pravdo & Shaklan (2009) announced a massive 
planet around the M-dwarf vBlO discovered by means of as- 
trometrical data. Zapatero Osorio et al. (2009; hereafter ZO09) 
made use of our NIRSPEC data set and found evidence for 
RV variations, which supported the planet hypothesis. They 
achieved a RV precision of about 300 m s However, this 
planet was later refuted by different groups: Bean et al. (2010b), 
who took high-resolution spectra (R = A/AA ~ 100,000) with 
CRIRES and who achieved a RV precision of ~ 10 m s ', and 
by Anglada-Escude et al. (2010). Additionally, Lazorenko et 
al. (2011) carried out an astrometric survey using the FORS2 
camera of the ESO/VLT on Cerro Paranal, Chile, but found no 
evidence for the existence of a massive planet orbiting vBlO. As 
part of this work, we aimed at finding out what had caused the 
spurious RV variations in the data analysis of ZO09. 

Here, we report relative RV measurements of 8 late M- 
dwarfs with NIRSPEC, and we support the capability of this 
instrument to detect giant planetary companions with short or- 
bital periods. In Section 2 we describe our M-dwarf sample, our 
observations and data reduction. In Section 3 we outline the de- 
tails of the data analysis, followed by the results and discussion 
(Section 4). 

2. Observations and Data Reduction 

As part of our M-dwarf survey (Deshpande et al., in prep.), 
we observed 8 M-dwarfs (2M2331, GJ1156, GJ406, GJ905, 
LHS1363, RXJ2208.2, and vBlO) at two or more epochs 
(Table [1]) using the NIRSPEC instrument, mounted on the Keck 

11 telescope on the summit of Mauna Kea in Hawaii (McLean 
et al. 1998). We aimed at conducting RV precision tests, and 
searching for RV drifts which could be interpreted as massive 
planets orbiting those M-dwarfs. Our sample comprised dwarfs 
with spectral types of M5.0-M8.0 and masses between 0.14 and 
0.075 Mq. Table |2] provides a list of the spectral types, /-band 
magnitude and the projected stellar rotational broadening vsin/ 
of the targets. 

NIRSPEC is a cross-dispersed, cryogenic echelle spectrom- 
eter employing a 1024 x 1024 ALADDIN InSb array detector 
In the echelle mode, we selected the NIRSPEC-3 (7-band) fil- 
ter and an entrance slit width of 0.432" (i.e. 3 pixels along the 
dispersion direction of the detector), except for the 2001 June 
observations of vBlO, where we used an entrance slit width 
of 0.576". The corresponding spectral resolving powers were 
R = A/AA ^ 22,700 and « ^ 17, 800, respectively for the 
0.432" slit and the 0.576" slit. The length of both slits was 
12". All observations were carried out at an echelle angle of 
~ 63°. This instrumental setup provided a wavelength coverage 
from 1.148 to 1.346 jum split into 10 different echelle orders, a 
nominal dispersion ranging from 0.164 (blue) to 0.191 A pix"' 
(red wavelengths). Weather conditions (seeing and atmospheric 
transparency) were fine during the observations, except for the 
2008 epoch, which was hampered by cirrus and strong wind. 



Table [T] lists the individual exposure times and the signal-to- 
noise ratios (SNRs) on average per spectral pixel in the stellar 
continua for each observing epoch. 

For each target, the spectra were collected at two different 
positions along the entrance slit. This nodding allowed later the 
removal of the OH sky emission lines. For the identification of 
atmospheric telluric absorption, near-infrared featureless stars of 
spectral types A0-A2 were observed close in time (on average 3 
min before or after the target observations) and position to our 
targets. 

Raw data were reduced using the echelle package within 
IRApQ. Nodded images were subtracted to remove sky back- 
ground and dark current. White light spectra obtained with the 
same instrumental configuration and for each target observation 
were used to flat-field the data. By adopting the apall task, 
we first identified and optimally centered the echelle orders in 
the two individual nodding frames for each target and traced 
these orders by adopting a second-order Legendre polynomial 
along the dispersion axis. In the next step, we extracted the one- 
dimensional spectra for each echelle order adopting the same 
aperture/trace parameters for both, the target as well as an arc- 
lamp exposure of Ar, Kr, and Xe, which was always acquired 
after observing the target and before pointing the telescope at 
the next star The air wavelengths of the arc lines were identified 
using the NIS10 database, and we produced preliminary wave- 
length calibration fits using a third-order Legendre polynomial 
along the dispersion axis and a second-order one perpendicular 
to it. The mean rms of the fits was 0.03 A, or 0.7 km s"'. 

3. Relative Radial Velocity Method 

We measured the RV of the stars relatively to the telluric lines 
present in the spectra as well as to a selected epoch of the star, 
employing a self-calibration approach. The radial velocity of the 
telluric lines is constant in all wavelengths down to a level of 
10 m S-' (e.g., Figueira et al. 2010, Seifahrt & Kiiufl 2008), 
which is about a magnitude smaller than the velocity precision 
we can achieve with NIRSPEC. The basics of the self-calibration 
method have been extensively described (e.g. Valenti et al. 1995, 
Endl et al. 2000, Bean et al. 2010a), so that we just give a con- 
cise description of the method here and point out the important 
aspects of its implementation. 

Briefly, the main idea is to model the observations and 
thereby determine the relative RV shift, and perform a fine tun- 
ing of the wavelength solution at the same time. Basically, the 
model spectrum is the product of a high-resolution telluric spec- 
trum with a Doppler-shifted version of a high-resolution refer- 
ence spectrum of the star. This product of those two spectra is 
then subjected to a convolution with the instrumental profile (IP) 
of the spectrograph, and finally binned to the sampling of the 
observed data. By variation of the free parameters of the model 
(Table O, the best fit model is evaluated by statistics. The 
input Doppler-shift which yields the best fit represents the mea- 
sured RV. 

Since our method requires the presence of telluric lines in 
the spectra, we restricted the analysis to the echelle orders 66, 
60, 58, and 57, which were heavily contaminated mainly by 
absorption lines of water vapor. These four orders correspond 

^ IRAF is distributed by the National Optical Astronomy Observatory 
(NOAO), which is operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National Science 
Foundation in the USA. 

^ http://physics.n ist.gov/PhysRefData/ASD/lines-form.htmll 
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Table 1. Journal of target observations. The average SNR values 
in the stellar continua are given. 



1.2710 1.2715 1.2720 1.2725 1.2730 1.2735 



Obs. date 


UT 


Sht 




Exp 


.(s) 


SNR" 


2MJ2331-2749: 




















2007-Iun-24 


14:33 


0.432 


X 


12 


2 


X 


200 




50 


2007-Iun-25 


14:12 


0.432 


X 


12 


2 


X 


200 




60 


GJ 406: 




















2007-Apr-30 


7:03 


0.432 


X 


12 


2 


X 


120 




50 


2007-Dec-23 


15:32 


0.432 


X 


12 


2 


X 


30 




50 


GJ 905: 




















2007-Iun-25 


14:42 


0.432 


X 


12 


2 


X 


20 




110 


2007-Oct-27 


10:51 


0.432 


X 


12 


2 


X 


120 




280 


GJ 1156: 




















2007-Iun-24 


7:24 


0.432 


X 


12 


2 


X 


120 




140 


2007-Jun-25 


15:55 


0.432 


X 


12 


4 


X 


300 




240 


LHS 1363: 




















2007-Oct-26 


12:09 


0.432 


X 


12 


2 


X 


300 




110 


2007-Oct-27 


12:13 


0.432 


X 


12 


2 


X 


300 




110 


LP 412-31: 




















2007-Oct-26 


12:44 


0.432 


X 


12 


2 


X 


300 




60 


2007-Oct-27 


12:56 


0.432 


X 


12 


2 


X 


300 




50 


RXJ2208.2: 




















2007-Iun-24 


13:49 


0.432 


X 


12 


2 


X 


120 




70 


2007-Iun-25 


13:40 


0.432 


X 


12 


2 


X 


120 




70 


vBlO: 




















2001-Iun-15 


14:06 


0.576 


X 


12 


2 


X 


100 




60 


2001 -Nov- 2 


4:43 


0.432 


X 


12 


2 


X 


120 




60 


2001 -Nov- 2 


5:39 


0.432 


X 


12 


2 


X 


120 




60 


2007-Jun-25 


13:22 


0.432 


X 


12 


2 


X 


120 




70 


2008-Jul-28 


6:07 


0.432 


X 


12 


2 


X 


120 




20 




" SNR on average in the pseudo stellar continuum per spectral pixel 
around 1.265 ^m. 



1.2710 1.2715 1.2720 1.2725 1.2730 1.2735 
wavelength [micron] 

Fig. 1. Comparison between the observed telluric spectrum 
(points), and the theoretical model (line). The observed telluric 
spectrum was taken by using a featureless A-star (HD181414), 
which was observed with NIRSPEC in the /-band on 2007-06- 
25. The rms of the telluric model fit to the observed data is about 
1%. 



Table 2. Properties of the M-dwarfs 



Name 


Sp. type 


J 


V sin / 

(km s-') 


Ref. 


2MJ233 1-2749 


M7.0 


11.65 


< 


12 


Desll 


GJ406 


M5.5 


7.09 




3 


ReilO 


GJ905 


M5.0 


6.88 


< 


3 


ReilO 


GJ1156 


M5.0 


8.52 


17.2 


+ 2.9 


Desll 


LHS 1363 


M6.5 


10.48 


< 


12 


Desll 


LP412-31 


M6.5 


10.48 


17.6 


±3.2 


Desll 


RXJ2208.2 


M5.0 


10.60 


18.6 


±2.3 


Desll 


vBlO 


M8.0 


9.91 


6 


5 


Moh03 



Abbreviations: 

Desll Deshpande et al., in prep. 

Moh03 Mohanty & Basry (2003). 

ReilO Reiners et al. (2010). 



to the wavelength ranges of A ~ 1.147 to 1.163 //m, 1.261 to 
1.279 //m, 1.304 to 1.323 fim, and 1.327 to 1.346 fim, respec- 
tively. For the echelle order numbering we refer to McLean et 
al. (2007). Subsequently, we subdivided each spectral order into 
5 equidistant pixel chunks of 200 pixel each (i.e. for all four or- 
ders together we have 20 chunks). This step was done to simplify 
the process of improving the model, to speed up the calculations 
and to account for variations of the IP throughout each spectral 
order. Each of the following steps was carried out on each chunk 
individually, and the SNR was determined by 

SNR = 5*/ a/^* + 5bg + Sbg2 + RON^ x 2n, (1) 



where 5* denotes the signal level from a star in electrons, in- 
tegrated over an aperture of n pixels, S bg is the signal level of 
the sky background, 5'bg2 is the signal level from the sky back- 
ground of the frame taken at the other nodding position, and 
RON denotes the read-out-noise level per pixel in rms electrons 
(for NIRSPEC, RON=65e '). The noise eiTors were propagated 
in the following data analysis steps. 

3.1. Step 1: Telluric template spectrum and determination of 
the instrumental profile 

For the calculation of the atmospheric transmission spectrum, 
we used the Line-By-Line Radiative Transfer Model (LBLRTM) 
code, which is based on the FASCODE algorithm (Cloiigh et 
al. 1992). LBLRTM is available as fortran source codeQ and 
runs on various platforms. As molecular database we adopted 
HITRAN (Rothman et al. 2005), which contains the 42 most 
prominent molecules and isotopes present in the atmosphere 
of the Earth. Following the approach presented by Seifahrt et 
al. (2010), we created a high-resolution theoretical telluric spec- 
trum for each observed spectrum by accounting for the air mass 
of the star as well as the weather conditions (water vapour den- 
sity column, temperature and pressure profiles) during the obser- 
vations. We retrieved the weather information from the Global 
Data Assimilation System (GDAS). GDAS models are available 
in 3 hours intervals for any location around the global. 



Source code and manuals are available under 
http : //rtweb . aer . com /lblrtm_description . html 
^ GDAS webpage: http : / /ready . arl . noaa . gov/READYamet . php 
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To calculate a first version of the instrumental profile (IP) of 
the spectrograph, we made use of the A-star observations next 
to our targets. First of all, we normalized the spectrum in such a 
way that the flux in the telluric continuum was at one. Next, we 
refined the wavelength solution of the observed telluric spectrum 
with the appropriate high-resolution theoretical telluric spectrum 
by adopting a second order polynomial. We then determined a 
preliminary version of the IP as the sum of 7 Gaussian profiles in 
a similar way as described in Valenti et al. (1995): Around a cen- 
tral Gaussian we grouped 3 Gaussians on each side of it, which 
allow to account for asymmetries in the IP. Free parameters were 
the height and width of the central Gaussian, plus the heights of 
the six satellite Gaussians (c.f. TableO. To reduce the number of 
free parameters per chunk, and to ensure that the method works 
robust, the positions and the widths of these satellites were fixed 
and set a priori in such a way that their half-widths overlapped. 

Next, we convolved the high-resolution theoretical spectrum 
with the determined preliminary IP and compared the resulting 
spectrum with the observed A-star spectrum. For a few telluric 
lines, we realized that the theoretical spectrum under- or over- 
estimated the line-depths. To produce a better match between 
theory and observation, we iteratively carried out a fine tuning 
of the line-depths in the high-resolution theoretical telluric spec- 
trum, then again convolved the modified telluric spectrum with 
the IP and evaluated the result with the observation by means of 
;\f^-statistics. The iterations were carried out until the reduced;^^^ 
reached 1. Fig. [T] shows a comparison between an A-star spec- 
trum (HD 181414) and the fit of the refined theoretical telluric 
spectrum as well as the residuals of the model fit to the A-star 
spectrum. The rms of the telluric model fit to the observed data is 
about 1% on average. The refined high-resolution telluric model 
spectrum from now on served the purpose of the telluric template 
spectrum. In the final step, we refined the wavelength solution of 
the observed spectrum, and then re-calculated the IP by adopting 
this new telluric template spectrum. 

3.2. Step 2: Stellar template spectrum 

Due to lack of appropriate theoretical model spectra which fit 
the stellar absorption features in the 7-band, we created the stel- 
lar template spectrum for one selected reference epoch by cal- 
culating an IP-free and telluric-free version of the target spec- 
trum. Concerning the reference epoch, we selected that epoch 
in which the stellar spectrum showed the highest SNR. To pro- 
duce the stellar template, we first applied the refined wavelength 
solution of the A-star spectrum (which was taken - on average 
- 3 minutes before or after the target observations) to the ob- 
served target spectrum of the same epoch. Since the telluric lines 
were present in the target spectrum, we needed to remove them 
from the spectrum. In preparation for this, we convolved the ap- 
propriate theoretical telluric spectrum with the IP. Then, we di- 
vided the target spectrum by the convolved theoretical telluric 
spectrum (Fig. [TJ. Similar to Bean et al. (2010a) and Blake et 
al. (2010), we found that this approach led to smaller uncertain- 
ties than when the usual method of the telluric lines removal was 
carried out, where the target spectrum is simply divided by the 
appropriate A-star spectrum. 

To create the final stellar template spectrum, we deconvolved 
the telluric-free target observation by the IP by employing the 
maximum-entropy method (MEM) with 5 times oversampling 
of the output spectrum. In the final step, we applied the refined 
wavelength solution that we had obtained for the A-star spec- 
trum to the 5-times over-sampled IP-free stellar spectrum, which 
from there served the purpose of the stellar template. 



1.2650 1.2655 1.2660 1.2665 1.2670 1.2675 1.2680 




1.2650 1.2655 1.2660 1.2665 1.2670 1.2675 1.2680 
wavelength [micron] 



Fig. 2. Example model components and fit for the radial velocity 
measurements. The components are given in the two top panels: 
the spectrum of the high-resolution theoretical telluric spectrum 
(top), and the deconvolved and RV-shifted version of the telluric 
free stellar spectrum (bottom). We note that the scale of the flux 
is different in each panel for better visibility. In the lower panel, 
we show the observed spectrum (points) and the best-fit model 
(line). 

Table 3. Free parameters in the model per chunk. 



Parameter 


degree of 




freedom 


Stellar absorption line depth 


1 


Linear stellar continuum trend 


2 


Doppler shift of stellar template 


1 


Telluric absorption line depth 


1 


Amplitude and width of main Gaussian 


2 


Amplitudes of satellite Gaussians 


6 


2'"' order wavelength solution vector 


3 



3.3. Step 3: Fitting the observed data 

For each target, we first determined the barycentric velocity dif- 
ferences Av'bc; for all observation epochs t with respect to that 
one of the stellar template epoch. This correction was calculated 
by use of the JPL ephemeris DE200 (Standish 1990). 

We constructed the model of the observation by multiplying 
the telluric template with the stellar template, where its Doppler 
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shift is one of the free parameters. The resuhing combination 
spectrum was subjected to a convolution with the IP that was 
determined in Step 1, and a new wavelength solution was cal- 
culated. Subsequently, all the free parameters (i.e. line-depths 
in the models, IP, See Table [3]l were refined by employ- 
ing Brent's optimization algorithm, and the fit to the observed 
data was evaluated by using statistics. The search range for 
the Doppler shift was Avbc,r ±15 km s"' with a step width of 
10 m s"'. We note that such a large interval would also allow us 
to detect large relative RV variations due to unseen massive com- 
panions such as low-mass stars and brown dwarfs. We calculated 
the values for each Doppler shift and then determined the 
exact ^^-minimum by using a Gaussian fit. That Doppler shift 
which led to the overall best fit model (^'^-minimum) constituted 
the measured RV of the star in the chunk, relatively to the stellar 
template. 

To determine the global (i.e. all chunks together) RV mea- 
surement, we combined all the RV measurements in all chunks 
into one by considering the following restrictions: Each chunks 
was given a specific weight which was determined from the 
average SNR in the stellar continuum, plus the number of tel- 
luric lines and stellar absorption lines which were present in 
that chunk, plus the depths of the stellar lines. Furthermore, 
we rejected chunks in which the RV-measurement constituted 
clearly an outlier (3cr above / below average of all RV measure- 
ments) by adopting sigma-clipping. No chunks were rejected 
for the stars GJ905, GJ1156, LHS1363, RXJ2208.2, and vBlO. 
For 2MJ233 1-2749 and LP412-31, one chunk was rejected each, 
while for GJ406 two chunks were rejected. All those rejected 
chunks were located in noisy areas with SNR levels lower than 
40 on average. We attribute these spurious RV shifts to improper 
stellar templates which contained artifacts coming from the de- 
convolution of low-SNR data. 

The global RV measurement was then determined as the 
arithmetic weighted mean of the un-rejected chunks. The error 
of the global RV measurement was determined as the weighted 
standard deviation of the un-rejected RV measurements in the 
chunks. 

4. Results and Discussion 

We analyzed the data sets with our relative radial velocity mea- 
surement approach and determined the relative RV measure- 
ments with respect to the selected reference epoch. For any of 
the eight M-dwarfs in our sample, we have not found significant 
evidence of relative RV variations at the level of 3cr (Table |4]i, 
where cr stands for the measurements uncertainty. The RV pre- 
cisions are in the order of 180-300 m s"', except for the obser- 
vations in July 2008, which were taken at low SNR. 

We investigated the period and mass range of companions 
which could be detected with such RV precisions. We deter- 
mined the minimum mass of the planet by employing a Monte- 
Carlo analysis, thereby probing planetary orbits with different 
parameters and investigating how many of these orbits could be 
recovered for the five measurements of vBlO. We considered 
only the case of a circular orbit and the mass of vBlO, which 
is OT* = 0.078 Mq. Fig.[3]shows the 3cr detection limit. We find 
that for companions with only a few days period, even planets 
with minimum masses of sin / > 0.3 Mjup can be detected 
with a RV precision of ~ 220 m s 

In Fig. 4 we show our relative RVs of vBlO and the mea- 
surements by ZO09. We note that for a proper comparison, we 
adopted the same reference epoch as in ZO09. The agreement 
between ZO09 and our measurements is within 1 cr of the quoted 
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D. 
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Fig. 3. Monte-Carlo analysis for the five vBlO measurements. 
For the mass of vBlO, we adopted 0.078 M© from Pravdo 
& Shaklan (2009). It is shown that with a RV precision of 
~ 220 m s"' even hot Jupiters with minimum masses Wp sin / > 
0.3 Mjup could be detected around late M-dwarfs with 3cr con- 
fidence. We note that for a larger number of measurements the 
number of aliasing peaks can be significantly decreased. 



uncertainties for all epochs except for the 2001 epoch (BJD = 
2452076). We provide next an explanation for the discrepancy of 
this one measurement. Similar to our data analysis, ZO09 used 
the telluric lines present in the target spectra as a stable refer- 
ence, but contrary to our analysis they did not account for any 
IP variations in their analysis, but calculated the RVs by cross 
correlation. In our analysis, we do not see any RV shift exceed- 
ing the RV-precision for any measurement. We get evidence that 
the diff'erent instrumental setting used on 2001-Jun-15 (0.576" 
slit instead of the standard setting of 0.432") produced an asym- 
metric instrumental profile (Fig. |5]l, which led to a significant 
RV shift when a simple cross correlation is adopted for the RV 
determination. Our results clearly demonstrate the importance 
of modeling the IP especially when observations are carried out 
with different instrumental settings. 

We compare our results to the work of Blake et al. (2010), 
who searched for companions to M- and L-dwarf wby using 
NIRSPEC at a spectral resolving power of ~ 25,000 in the K- 
band. They adopted one spectral order covering the wavelength 
range from 2.285 to 2.318 [im. to measure the dense and strong 
CO-absorption line pattern present in those dwarfs. As a stable 
wavelength reference, they made use of the CH4 telluric absorp- 
tion lines present in the observations. Similarly to us, they em- 
ployed a self calibrating approach, with the diff'erence that they 
adopted theoretical models for M- and L-dwarfs, which well- 
described the observations. 

Blake et al. obtained measurements with SNRs in the range 
of 50 to 100 in the pseudo stellar continua, and they report 
RV precisions of 100-300 m/s for slowly rotating late-M and 
L dwarfs. The uncertainty of Blake et al. in the K-band is in 
agreement with our derivation of 180-300 m/s in view of our 
SNRs. However, our wavelength coverage is about twice that in 
Blake et al. According to the relative RV precision formulae, we 
should have obtained better velocity precision in terms of wave- 
length coverage, which is not the case. We conclude that both 
the larger number of deep lines (more than 30 lines with a line 
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Table 4. Relative radial velocity measurements. 
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Fig. 4. Our RV measurements of vBlO (crosses) vs. the measure- 
ments of ZO09 (open circles). The RV precision given by our 
analysis is about 220 m s \ except for the last epoch in 2008, 
which was hampered by bad weather. We prove that the RV shift 
in the work of ZO09 in the first 2001 epoch (BJD= 245 2076) 
originates from unaccounted asymmetries in the IP rather than 
from a planetary companion. 
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Fig. 5. Instrumental profiles (IPs) of NIRSPEC for two different 
observing epochs of vBlO. On 2001-06-15, a broader slit was 

used (0.576"; solid line) than for the reference epoch (0.432"; 
3007-06-25; dashed line). To visuahze the asymmetries between 
both IPs, we calculated the ratio between both IPs and scaled 
the resulting function for better visibility (dotted line). ZO09 did 
not account for these asymmetries between both IPs, which led 
to a spurious RV shift of about 1 km s"^ for the 2001-06-15 
measurement in their data analysis. 



depth larger than 50%) in the CO-band region as compared to the 
J-band (only a few lines with a depth larger than 25%) as well 
as the use of theoretical template spectra instead of deconvolved 
stellar spectra appear to compensate for the shorter wavelength 
coverage in a similar factor (c.f. Equation 6 in Butler et al. 1996). 

We note that Reiners et al. (2010) and Rodler et al. (2011) 
carried out theoretical RV precision studies of M- and L-dwarfs, 
by adopting theoretical models of M-dwarfs (e.g. del Burgo et 
al. 2009). As result, they find that the highest RV precision for 
M-dwarfs is attained in the Y band around 1 yum, rather than 
in the J- , H- or /T-band. For L-dwarfs, however, Rodler et 



Obs. Date 


BJD 


rel. RV 




245 0000-1- 


(m s-i) 


2MJ2331-2749 






2007-Jun-24 


4276. 10874 


194 ± 224 


2007-Jun-25 


4277.09424 


ref. epoch 


GJ 406: 






2007-Apr-30 


4220.79751 


reference epoch 


2007-Dec-23 


4458.14962 


-11 + 238 


GJ 905: 






2007-Jun-25 


4277.11249 


-233 ± 201 


2007-Oct-27 


4400.95657 


reference epoch 


GJ 1156: 






2007-Apr-30 


4220.81513 


141 ± 185 


2007-Dec-22 


4457.16844 


reference epoch 


LHS 1363: 






2007-Oct-26 


4400.01234 


-11 ± 196 


2007-Oct-27 


4401.01511 


reference epoch 


LP 412-31: 






2007-Oct-26 


4400.03654 


reference epoch 


2007-Oct-27 


4401.04491 


298 ± 260 


RXJ2208.2: 






2007-Jun-24 


4276.07879 


-189 + 307 


2007-Jun-25 


4277.07266 


reference epoch 


vBlO: 






2001-Jun-15 


2076.08951 


19 + 230 


2001-NOV-02 


2215.70560 


69 ± 223 


2OOI-N0V-O2 


2215.74225 


-3 ±217 


2007-Jun-25 


4277.05865 


reference epoch 


2008-Jul-28 


4675.75669 


131 ±497 



al. (2011) reported that the highest RV precision is attained in 
the 7-band. 

We conclude that for an accurate relative RV determination 
with NIRSPEC, a self-calibrating approach, which accounts for 
changes in the instrumental setting, produces the best measure- 
ments in terms of RV precision. Although with our RV precision 
we would be able to detect massive hot Neptunes around late 
M-dwarfs, we have not found any brown dwarf or massive plan- 
etary companion in our survey. Additionally, the re-analysis of 
the data of the M8-dwarf vBlO presented in ZO09 now clearly 
confirms the non-existence of a massive planet orbiting that 
dwarf and agrees with the results by other research groups (e.g. 
Anglada-Escude et al. 2010; Bean et al. 2010b; Lazorenko et 
al. 2011). 
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